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Abstract 

CN ■ 
O ■ 

' We analytically derive transition probabilities for four-neutrino oscillations 

■ in matter. The time evolution operator giving the neutrino oscillations is ex- 
Q . pressed by the finite sum of terms up to the third power of the Hamiltonian in 

\ a matrix form, using the Cayley-Hamilton theorem. The result of computa- 

Q>^ \ tion for the probabilities in some mass patterns tells us that it is realistically 

difficult to observe the resonance between one of three active neutrinos and 

■ the fourth (sterile) neutrino near the earth, even if the fourth neutrino exists. 

i> : 

a^ ■ 

O ; I. INTRODUCTION 

o . 

■ A neutrino oscillation is a transition among neutrino flavors. Several types of the ob- 
servations tell us that neutrino oscillations occur [1-8]. They are classified into the solar, 
atmospheric and LSND experiments. 

^ , The mass squared differences are the parameters showing the neutrino oscillations. In 

g_) \ order to explain three kinds of neutrino experiments within one framework, three kinds of 

■ the mass squared differences are needed. Therefore we consider the four-neutrino oscillation, 
' where the fourth neutrino doesn't have the weak interaction. Three active neutrino flavors 

^ . (i/e, t'/i, i^t) interact with leptons in the weak interaction. So the fourth neutrino is called 

^ ! sterile neutrino (i/^). 

" " ' The neutrino oscillation pattern in vacuum can get modified, when the neutrinos pass 

through matter. It is known as the Mikheyev-Smirnov-Wolfenstein (MSW) effect [9], which 
can be described by an effective Hamiltonian. The interaction with the neutral currents 
occurs for three active neutrinos. Thus, for the three active neutrinos, one doesn't need to 
consider the interaction with the neutral currents [10]. But the sterile neutrino has neither 
the charged- nor the neutral-current interactions. It means that one needs to consider the 
effect of the matter interacting with the sterile neutrino and to introduce the 4x4 mixing 
matrix of four neutrinos which is an extension of the 3x3 Maki-Nakagawa-Sakata (MNS) 
matrix [11]. 

Analytical calculations of active three-neutrino oscillations in matter have been per- 
formed [12]. In this article, we derive analytically the transition probabilities for four- 
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neutrino oscillations. Our calculations include the effects of the interaction with charged 

and neutral currents. 

The outline of the article is as follows. In Sec. II, two kinds of the bases to express 
four-neutrino states are introduced. These bases are connected by a mixing matrix. To 
describe neutrino oscillations in matter, the effective Hamiltonian with charged and neutral 
current is introduced. In Sec. Ill, we calculate the transition probabilities from the effective 
Hamiltonian. In order to derive the transition probabilities, we make use of the Cayley- 
Hamilton theorem and the formula for the root of biquadratic equation. In Sec. IV, the 
transition probabilities are concretely computed in two cases of the four-neutrino oscillation 
schemes. Finally, in Sec. V, we discuss the effects of four-neutrino oscillations in matter. 



II. FORMALISM 

A. Two bases and a mixing matrix 

Neutrinos are produced in flavor eigenstates \i'a){oi = e, /x, r, s). Between the source and 
the detector, the neutrinos evolve as mass eigenstates |i^a)(a = 1,2,3,4). There are two 
kinds of eigenstates: and jz^a). These eigenstates are defined by neutrino fields Va ^-nd 
z/q corresponding to each eigenstate: I'^O) = |z/), |z/q) = \a), l^a) = |a), where a vacuum 
state is given by |0). In the present analysis, we will use the plane wave approximation of the 
fields. In this approximation, a neutrino flavor fleld Uq. is expressed by a linear combination 
of neutrino mass field Ua'. 

4 
o=l 

where C/ is a 4 x 4 unitary matrix with the elements Uaa- If we write this relation in neutrino 
eigenstates, then 

4 
a=l 

An arbitrary neutrino state ip is expressed in both the flavor and mass bases as 

4 4 / \ 4 

a=e,ix,T,s a=e,ix,T,s a=l a=l \a=e,/i,T,s / a=l 

where ■^^ and are the components of of the flavor eigenstate basis and the mass 
eigenstate basis, respectively. They are related to each other in the form with respect to 

^a= Y (4) 

a=e,ix,T,s 

If we deflne the matrix elements as 
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U^l Ufj,2 Uij,3 
UtI Ut2 Ut3 UrA 
\ Usl Us2 U,3 UsA J 



^3 



(5) 
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the relation between the flavor and the mass eigenstates are 
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The unitary matrix U is the mixing matrix of four neutrinos. There are 6 mixing angles and 
3 mixing phases as parameters of C/, in the case of 4 neutrinos. In this analysis, we ignore 
the CP violation by putting the mixing phases equal to zero. Then, U is the real orthogonal 

matrix [12]. 

A parameterization iov U — U {9i2, ^13, ^14, ^23, ^24, ^34) is given by 



U = C/34C/24C/14C/23C/13C/12, 



(8) 



where matrix elements are 



Cij — COS 6ij , Sij — sin 6ij , 



(9) 
(10) 



and the mixing angles 6*12, 6*13, 6*14, 6*23, 6*24, 6*34 [13,14]. By this definition, the mixing matrix 
becomes 
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u = 



CiaS 



14^13 



'14 



-C23C24SU 

-C12C24S13S23 

-C12C13S14S24 



C12C23C24: 

— C24S12S13S23 

— C13S12S14S24 



C13C24S23 
— S13S14S24 



C14S24 



—C12C23C34S13 

+C34S12S23 

—C12C13C24S14S34 

+C23S12S24S34 

13'S'23'S'24'S'34 



— C23C34>S'i25'i3 
—C12C34S23 
—C13C24S12S14S34 
—C12C23S24S34 

13'S'23'S'24'S'34 



Cl3C23C'34 

— C245'i3S'i4S'34 

— Ci35'235'24<S'34 



Ci4C24S'34 



— Cl2Cl3C24C'34'S'l4 — Cl3C24C'34>S'l2>S' 



a2L'13<-^24<-^ 34^14 
+ C23 C34 <S'i2 <S'24 
+ (^12(^34 'S'135'23 5*24 
+C*12C'23'S'l35'34 
— 5'i2<S'23<S'34 



13l_/24<-^34^12^14 
— C12 C23 C34 S 24 
+C34'S'i25'i3S'235'24 
+ C'23'S'l2'S'i35'34 
+Ci2<S'23<S'34 



-C24C'345'l3S'i4 

-C13C34 5*23 5*24 C'l4C'24C'34 

-Ci3C235'34 



B. Hamiltonian in matter 



In the mass eigcnstatc basis, the Hamiltonian TIq participating in the propagation of 
neutrinos in vacuum is given by 



f El \ 

^2 

^3 

\ E4 / 



(12) 



where Ea {a = 1, 2, 3, 4) are the energies of the neutrino mass eigenstates \a) with mass ma- 

Ea = y/rria^+P''- (13) 

Here and hereafter we assume the momentum p to be the same for all mass eigenstates. 

There are two kinds of the additional potentials for describing the interactions between 
neutrinos and matter. One is the interaction of the charged particles (electrons) and its 
neutrino Ve'- 

Ve = V2Gf diag(7Ve, 0, 0, 0). (14) 
The other is the interaction of the neutral particles(e.g. the neutron) and active neutrinos 

(i/e,i^/i,^'T): 

Ko = V2Gf diag(-^7V„, -^7V„, -^7V„, 0), (15) 
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where Gp, Ng and Nn are the Fermi weak coupling constant, the electron number density and 
the neutral particle number density, respectively. Note that we assume the particle number 
densities to be constant throughout the matter where the neutrinos are propagating. 
The interaction term (15) can be separated into two parts as 



Vn^V2GF diag(0,0,0,+^7V„), 

V = V2Gf diag(-^Ar„, -^Nn, —Nr., —N^). 



(16) 
(17) 

(18) 



These interaction terms are written by the flavor eigenstate basis. Therefore the interaction 
terms in the flavor eigenstate basis must be transformed into those in the mass eigenstate 
basis by the mixing matrix U. The interaction terms in the mass eigenstate basis are 



U-^VM = A 



U-^V„U = A. 



Ue2Uel ^62^43 ^^62^^64 

\ UeiUel Ue4:Ue2 Ue^UeZ U'^A 

( UslUs2 UslUsS UslUsA \ 

Us2Us\ Us2Us3 Us2UsA 

\ UsiUsl UsiUs2 UsaUsZ 



U'Wu -- 

where / is the 4x4 unit matrix, and the matter densities A^, An and A are deflned by 



V2GfN^ 
1 



V2 



GpN„, 



A, 
1 



-A— 
2 



Thus, the Hamiltonian in the case when the neutrinos propagate in matter is 

nm^no + u-'vM + u-'VnU - aj. 



(19) 

(20) 
(21) 

(22) 
(23) 

(24) 



III. CALCULATIONS OF THE NEUTRINO TRANSITION PROBABILITIES 

The transition probabilities are represented by the time evolution operator. In the flavor 
state basis, the unitary transformation from the initial state 'ipf{t = 0) to the final state 
ipf{t) is given by the operator 

Uf{t) = Uf{t,0), (25) 

where [7/(^2, ^i) is the time evolution operator from time ti to t2 in the flavor state basis. 
The Hamiltonian Tiflavor in the flavor state basis is represented by using the mixing matrix 
U and the Hamiltonian Tim in the mass state basis: 



TCtla 



UHmU- 



(26) 
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The Schrodinger equation in the mass eigenstate basis is 

^^V'm(^) = nm^mit). (27) 

Equation (27) has a solution 

V'm(i) = e-^^-Vm(0), (28) 

where e"'^""* is the time evolution operator. Inserting t — L into (28), the solution of the 
Schrodinger equation (27) is 

iPUL) = ^m{t) = e-^^-^Vm(O) = Um{L)^m{0), (29) 

t=L 

where L stands for the distance through which neutrinos run for the time t, because the 
speed of neutrinos is almost equal to that of light. 

The neutrino state ■0/(-f') at i = L in the flavor state basis is expressed as 

^;(L) = UiJmiL) = t/e-'^'»^^^(0) = ^7e-^^-^t/-V/(0) = Uf{L)ijf{0). (30) 



A. Traceless matrix T 

In order to find the explicit form of the time evolution operator e"*^"**, which is the 
exponential of the matrix, the Hamiltonian in the matrix form is separated into the diagonal 
and the traceless matrices. The trace of the matrix Tim in (24) is 

trHm ^ Ei + E2 + E3 + E^ + Ae- SAn, (31) 

where we use the unitarity conditions, e.g. Ud'^ + Ue2^ + Ues^ + Uei' = 1. 
An arbitrary N x N matrix M can always be written as 

M = Mo + ^(trM)7jv, (32) 

where Mq and In are N x N traceless and unit matrices, respectively. Note that trMo = 0. 
Then the 4x4 matrix can be written as 

Hm^T+^itiUm)!, (33) 
where I4 — I and the matrix T is traceless. The matrix T can be written by 
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\ 

^2,134 

£;3,i24 
V £;4,i23/ 

(Uj-\ [/el[/e2 f/elf/.3 f/elf/e4 \ 
Ue2Uel Ue2 — J t^e2^^e3 ^^62^^64 



^e3 

UesUel Ue3Ue2 f^e3 ~ 4 f^e3^e4 



\ UeiUel Ue4Ue2 UeiUeS, Uei' — \ ) 

(UsX-\ UslUs2 UslUsS UslUs^ \ 

Us2Usl Us2 — 4 Us2Us3 Us2Us4 

UssUsl Us3Us2 UsZ — \ Us?JJs4: 

\ UsaUsI UsdJs2 Us4:Us3 Ugi ~ \ j 



(34) 



where E^b (a, 6 = 1, 2, 3, 4, a 7^ h) and Ek^imn{k, 1,171,11 — 1, 2, 3, 4) are defined as 

-£■06 = Ea — Eh, Ek^lmn = E^l + E^m + Ekni 

respectively. The energy differences Egj, are not linearly independent, since they obey the 
following relations: 



E. 



ab 



— Eba, E12 + £"23 + -£'31 = 0, £'12 + -E24 + -£'41 = 0, Ei3 + £34 + E. 



41 



0. 



Thus, only three of the EabS are linearly independent. 

Therefore the time evolution operator e"*'^"'* can be rewritten by the traceless matrix T: 



where 



UmiL) = e 
g-i(trWm)L/4 jg phase factor. 



—iHrnL 



(f)e 



-iTL 



(35) 



B. The Cayley-Hamilton theorem 

In order to find the concrete form of the definite matrix e~*^^, we use the Cayley- 
Hamilton theorem. The exponential of the 4x4 matrix T can be expressed by an infinite 
series 

e"^^^ = kol + hT + k2T'' + + kj"^ + • • • , 

where kn = ^; (n = 1, 2, 3, 4, • • • ). The Cayley-Hamilton theorem implies that the 
eigenvalue A in the characteristic equation 

det(T - A74) = + caA^ + C2A^ + ciA + cq = (36) 

of the matrix T can be replaced with T to give 

+ caT^^ + C2T^ + ciT + cq = 0, (37) 



7 



where Cj {j — 0, 1, 2, 3) are coefficients. Using (37) repeatedly, the matrix e can be 
formally written in the form 



-iTL 



(38) 



where aj {j = 0, 1, 2, 3) are coefficients which differ from A;„ and Cj in general. 

Because T is the definite matrix, we need to find the coefficients aj explicitly in order 
to obtain the matrix e~'^^^. If the characteristic equation (36) has four solutions {k — 
1, 2, 3, 4), one can write the eigenvalues of e~'^^ as 



e '^'^^ = ao + aiAfc + asA^^ + asA^l 



By defining the following matrices. 



-,—iL\2 
-,—iL\s 



A 



(\ -iLXi -L^Ai^ +iL^Xi^\ 
1 -iLX2 -L^X2^ +iL^X2^ 
1 -iLXs -L^Xs^ +iL^Xs^ 



a 



ai 

0.2 



\ J \ 1 -iLA4 -L'^X4^ +iL^X/ J 

(39) is written as the matrix form e = Aa. Then, one obtain the coefficient a by 

a = A"^e. 

Therefore, we should find the eigenvalues A^ of the matrix T in order to know a. 



(39) 



(40) 



(41) 



C. Characteristic equation of the matrix T 

In order to obtain the eigenvalues A^ of the matrix T, one must solve the characteristic 
equation 



Til -A 


T12 


Tl3 


Tu 


T21 


T22 — X 


^23 


T24 


T31 


T32 


T33 — A 


T34 


T41 


T42 


^43 


T44 — A 



A^ + C3A^ + C2A2 + ciA + Co = 0. 



(42) 



The coefficient Cq is the determinant of T, ci and C2 are expressed by the sum of cofactors 
of T, and C3 is given by the trace of T, 

C3 = -(Til + T22 + T33 + T44) = -trT = 0, 



C2 — -\-T11T22 + T'iiT33 + T11T44 + T22T33 + T22T44 + T33T44 

Ti rri rri rri rri rri rri rri rri rri rri rri 

12-^21 — -t 13-^31 — -t 14-^41 — -t23-'32 " -'24-' 42 " -'34-' 43; 
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T' rri rri rri rri rri rri rri rri rri rri rri 

ll-'22-'33 ~ J- 11 -1- 22-1-44: — -tll-'33-'44 — -'22-' 33-' 44 

+^11^23^32 + ^11^24^42 + 711734743 + T22Ti^T2i + T22T14T41 + T22Tz4T4-i 

+733712^21 + 733T14T41 + T33724T42 + T44T12T21 + T44Ti3T'3i + T44723732 

Ti rri rri rri rri rri rri rri rri rri rri rri 

13-t34-tll — -fl-l-t-i;U31 — -Ll2-L21-L'll — -tll-tl2-t21 

TI rri rri rri rri rri rri rri rri rri rri rri 

12-'23-'31 ~ -'13-' 32-^21 ~ -'23-' 34-^42 ~ -1 24-^ 43-^ 32, 



Co — 711T22733T44 

T' 'T"' T"' T"' T"' 

ll-'22-' 34-^43 ~ -f ll-f 33-* 24-1 42 " ll-f 44-4 23-^ 32 

Tr rri rri rri rri rri rri rri rri rri rri rri 

22-L 1 l-t 11 ~" ^ 22-L 1 l-t 31 ^ -t 33-^ 1 1-t 12-t 21 

+7^11723734742 + 7^11724743732 + 7227^13734741 + 7227^14743731 

+7337^12724741 + 7337^14742721 + 744T12T23T31 + 744T13T32T21 

TI rri rri rri rri rri rri rri rri rri rri rri 

12-' 23-' 34-' 41 ~ -'12-' 24-' 43-^31 ~ 13-* 32-* 24-* 41 

TI rri rri rri rri rri rri rri rri rri rri rri 

IS-t 34-1 42-1 21 — -t 14-1 42-1 23-^31 ~ -t U-t 43-^ 32-1 21 

+7l272l734743 + 7\373l7'24742 + 7i474l723732- 

The four roots of the biquadratic equation (42) are given from the solutions of the two 
quadratic equations [15] 



X^±Vh^2X + ^j + ^^j- -co = 0, (43) 

where to is one of real roots of the cubic equation 

- 02^ - Acot + 4coC2 - ci^ = 0. (44) 
Note that C3 = — trT = due to the definition of T. 

D. Calculation of time evolution operator 

Prom (35) and (38), the time evolution operator is written as 

= [aol + {-iLT)ai - L'^T^a2 + iL^T^as] . (45) 
The coefficients a^ 's are given by (41) as follows, 
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= -^2A3A4 ^-iLAi _ A1A3A4 ^-iLAa _ A1A2A4 ^-iLAg _ A1A2A3 ^_^2.A4 
A12A13A14 A21A23A24 A31A32A34 A41A42A43 

i I A2A3 + A2A4 + A3A4 _iLX^ . A1A3 + A1A4 + A3A4 _iLX-> 
1j \ A12A13A14 A21A23A24 

^AiA2 + A1A4 + A2A4^_^^^g ^ A1A2 + A1A3 + A2A3^_^^^^ 

A31A32A34 A41A42A43 

J_ f A2 + A3 + A4 ^_^^^^ _^ ^]_+^3_+A4^_iLA2 

\ A12A13A14 A21A23A24 

^ Ai + A2 + X4 ^-iLX3 Ai + A2 + ><3 ^-iLX4 
A31A32A34 A41A42A43 



°3 = 73 , , , g + X X X e + , , , e + , , , e- . , 

L, \A12A13A14 A21A23A24 A31A32A34 A41A42A43 / 

where Xab = Aa — A^. Inserting these results into (45), one can find the time evolution 
operator in terms of Aa's. For example, the term containing e"*^^^ is 

A2A3A4/ - (A2A3 + A2A4 + A3A4)T + (A2 + A3 + A4)r' - ^_iLx, ^4g^| 
Ai^ — (A2 + A3 + A4)Ai^ — (A2A3 + A2A4 + A3A4)Ai — A2A3A4 

Using the relations of the coefficients and solutions for the biquadratic equations, 

Ai + A2 + A3 + A4 = -C3 = 0, (47) 

A1A2 + A2A3 + A3A4 + A4A1 + A1A3 + A2A4 = C2, (48) 

A1A2A3 + A2A3A4 + A1A2A4 + A1A3A4 = -ci, (49) 

A1A2A3A4 = Co, (50) 

(46) can be written as 

(ci + C2A1 + Ai=^)7 + (C2 + Ai^)r + AiT^ + T^ ^_,Lx, 
4Ai^ + ci + 2C2A1 

Therefore, the matrix e"^^^ is given by 



(51) 



J]S„e-^^^», (52) 
0=1 

(ci + C2A„ + A „^)7 + (c2 + A„^)r + A„r2 + 
4A„^ + ci + 2c2A„ 



-"a = , X 3 , _ , r, X ■ W-^J 



The time evolution operator in the mass eigenstate is derived by 
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Um{L) = e-^^-^ = J]S„e-^^^». (54) 

a=l 

Using the mixing matrix [/, the time evolution operator in the fiavor eigenstate is given as 
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Uf{L) = e-''^f^ = Ue-''^"^^U-^ = <PY1 Bae~'^^% (55) 

. UB^U-^ = (c.+c.A„ + V)J+fe + V)f + A„f- + f3 

4A/ + C1+2C2A. ^ ' 

where T = UTU~^. 



E. Transition probabilities in matter 

A probability amplitude is defined as 

Aaf} = {/3\Uf{L)\a), a, /3 = e, /X, T, s. (57) 
Inserting (55) into (57) the probability amplitude becomes 

4 

A^p^<Pj2iBa)^pe-''^''^, (58) 

a=l 

. ^ (Cl + C2Xa + Xj)5af, + (C2 + Xa')f^fS + Xa\f\^ + {f^)^p 

^ "^"^ 4A„=^ + ci + 2c2A„ ' ^ ^ 

where 

= {a\m^f^p, (a|T2|/3) = (TV (a|r3|/3) = (T^^. 

Here 5^^, Tq,^, (T^)^^ and (T''^)^,^ are all symmetric. The probability of transition from 
the neutrino flavor a to the neutrino flavor (5 is defined by 

P^p = \A^p\^ = A^p*A^p. (60) 

Using the definition of the probability amplitude (57), one finds 

^a/^ = E E mapiB,)^,e'^'^^'-''\ (61) 

a=l 6=1 

where the symmetry of T leads P^p — Ppa- 

The probabilities for the oscillations in vacuum are given by setting A^ — A„ = in the 
definition of the probability amplitude. From (24), one can find 

{b\UmiL)\a) = (&|e-^^-^|a) = {b\e-''^°''\a) = e-'^'^''5ab, a, 6= 1,2,3,4, 

where Hml a _n a -n — ^0- Setting the following parameters 

Xab = -EabL, 6afi = E ^aa^^lBa^ + ^ E ^ocaU fiaUabU jib-, (62) 

a a<6 
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the probability amplitude A^p and the probability Pa/s of transition are given by 

4 

Aap = mUfiL)\a) = (/3|C/e-^«-^C/-» = J] C/aaC//3ae-'^»^, (63) 

a=l 

-Pa/3 = - 4 ^ U ^aU j3aU ahU ph sin^ (64) 



a<6 

The expression (55) gives the relations at L = 

4 



a=l o=l 



Then, we can rewrite the probabilities (61) for oscillations in matter in a form analogous to 
(64) in the vacuum, 



4 4 



PaP = - 4^^ Y^{Ba)^p{Bb)^p Sin^ Xab, (65) 
a=l 6=1 

where ^ 

Xah = ^(Aa - Aft). (66) 

If one set q; 7^ /9, then it is written as 

4 4 

o=l 6=1 
a<6 

2^ 



_ (^2 + A„^)r.;, + XgjTXp + (J^^).;3 

4Aa + Ci + 2C2Aa 

Prom unitarity, (65) gives the relations 



Pee + Pe^. + Per + Pes = 1, (69) 

P^e + Pt.^ + Pt.r + P^.s = l, (70) 

Pre + PrM + ^- + Prs = 1, (71) 

Pse + Ps,. + + Pss = l, (72) 

where 

-fe// -Pjue; Per Prei PfiT Pt/xi Pes Psei Pfis Psfxi Pts Pst- (73) 

Hence, there are only 6 independent transition probabilities. 
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IV. FOUR-NEUTRINO OSCILLATIONS IN MATTER 



We apply some results obtained in the previous section to the four-neutrino oscillation 
models. The probability (65) for transition of four-neutrino oscillation in matter contains 
the neutrino energy differences Eab — Ea — Ef, {a,h — 1,2,3,4) in vacuum. This energy 
differences are approximately given by the matter mass differences Am^^ = uia^ — rrib^, 
which are well-known quantities in various kinds of the neutrino experiments. From the 
on-shell condition (13), 

Am„;.2 = - m,^ = EJ - E,^ = {E^ - E,){Ea + E,) = 2E,, ^"^^\ (74) 

Now, we assume the average E of the neutrino energies is about lOGeV : 

E = Er + E, + E, + E, ^ ^^^^^ 
4 

One defines the difference S between two kinds of the averages of two neutrino energies, e.g., 
the average of the neutrino #1 and #2 and the average -Si±^ of the neutrino #3 and 

#4. We assume that the difference 5 = ^2±Ea _ Mi±M2. [g about leV which is estimated from 
the result of the LSND experiment. Then, 

- ) + 2- 



It is found ioT 6 <^ E that 

In four-neutrino oscillation analysis, there are three kinds of the mass squared differences. 
They are used as the parameters in the solar and atmospheric oscillations and the LSND 
experiment, which are represented as Amg^ig^j,, Aml^^ and AmLSNp, respectively. Using 
these mass squared differences, one can consider several distinct types of mass patterns [10]. 
They are classified into the so-called (3-|-l)-scheme and (2-|-2)-scheme. We concentrate the 
discussion on two of the several mass patterns in Fig. 1. The phenomenology and the mixing 
matrix depend on the type of the mass schemes. 



A. (3-|-l)-scheme 

We assume the mass patterns shown in Fig. 1(a), in which there are three close masses 
and one distinct mass. Let 777.4 and Aml^ be the distinct mass and the largest mass squared 
difference, respectively. 

First, three kinds of the neutrino mass squared difference are put as follows [12,14], 

Ami, = AmLa. ^ 10-^eV^ (77) 
Ami, = Am^tm ^ 10-'eV^ (78) 
Ami, = Am^sND - leV^ (79) 
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Then, the energy differences -Bab's are expressed by 




771 IP \ TP IP IP IP IP TP IP 

^31 — ^32 -r — — -C/43 — — -C/31, 



(80) 



where we suppose that the average E of the neutrino energies is lOGeV. 

Next, we consider the approximate mixing matrix for the (3+l)-scheme [14]: 




(81) 



where e is small: < e < 0.1. Here the 3x3 sub-matrix that describes the mixing of the 
three active neutrinos has the bimaximal form. The mixing matrix of (81) is given from (11) 
by taking 



As an illustration of the resonance phenomena, the energy differences |Aa — A^l (a, 6 = 
1, 2, 3, 4, a 7^ 6) are plotted as a function of the matter density A in Fig. 2. Here we assume 
that the electron number density A^e is eqiial to the neutral particle number density : 

= Nn. That is, = A and A„ = |A in (22) and (23). Note that Aa's are the effective 
neutrino energies in matter. In Fig. 2, the resonances occur when the energy levels in the 
presence of matter approach the values of each other. 

The transition probabilities (65) for the neutrino oscillations in matter as functions of 
the matter density A are shown in Fig. 3. They are some remarkable results showing the 
effects of the sterile neutrino. Here, we set a parameter 77 = 1. The parameter 77 is defined 
&s L/E = rj X (2i?/10GeV) = rj x 6.46 x lO^eV"^, where 2R is the diameter of the earth : 
R = 3.23 X lO^^gV-i = 6378km. 

From Fig. 3, the following results can be derived. If there is a httle mixing of the 
neutrino #1 and #4 in the neutrino mass states, i.e., e 7^ 0, the fourth (sterile) neutrino 
effects appear as resonance. For example, the probability Pee has a little transition beyond 
the matter density A ~ 10~^^eV and a sharp drop at ^4 ~ 10~^°eV. 



The mass pattern about the (2-|-2)-scheme is shown by Fig. 1(b). We assume that both 
mass differences of the Ue and Vg and of the and Vr are small. The three kinds of the 
neutrino mass squared differences are put as follows. 




(82) 



B. (2-|-2)-scheme 



Am^i = Amlsj^^ ~ leV^ 
Aml^ = Aml^ ~ 10-^eV^ 




(83) 
(84) 
(85) 
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where the average E of the neutrino energies is treated as lOGeV. 
The approximate mixing matrix for the (2+2)-scheme [14] is 
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(86) 



where e is supposed to be smaU: < e < 0.1. These parameters resemble those in Ref. [14], 
except for ^13 = 0. The mixing matrix in (86) for the (2+2)-scheme is given from (11) by 
taking 

Ou^e, ^13-0, ^24-0, ^34 = 0. (87) 



6^23 - ^, 



We take = 1, as shown in the (3+1) scheme. 

In the (2+2)-scheme, the results of the energy differences |Aa — \}\ {a,b = 1,2,3,4) are 
presented in Fig. 4 as a function of the matter density A. The transition probabilities Pap 
for the neutrino oscillations in matter are shown in Fig. 5. 

We assume that two mixings between and Ug and between i/^ and i/r are maximal, and 
the other four mixings are minimal as (87). Some transition probabilities Pap are shown in 
Fig. 5, and Pgr resembles Pg/^ in a probability pattern. The transition between two neutrinos, 
of which masses are clearly distinguished each other, occurs beyond A ^ 5 x 10~^^eV. 



V. DISCUSSION 

The main result of our analysis is given by the time evolution operator (55) for the four- 
neutrinos in matter. The time evolution operator (55) in the flavor eigenstate is expressed 
as a finite sum of elementary functions in the matrix elements of the Hamiltonian (24). The 
transition probabilities in matter have been given by (65). We also have analyzed the matter 
effects of the transition probabilities, assuming that there are four kinds of neutrinos and 
that the fourth (sterile) neutrino has a little mixing with the other neutrinos. 

The resonance between an active neutrino (ve, i^/i or Pr) and the sterile neutrino Ug occurs 
at the matter density A ~ 10~^''eV. Is this matter density realistic? We consider about the 
density of the sun (see Appendix A). Using a solar model [17], the electron matter density A 
in the sun is 1.06 x 10~^^eV at the center of the sun and 2.78 x lO^^^eV at the surface of the 
sun, respectively. The average of the electron matter density in the sun is 1.40 x lO^^^eV. 
Thus, the electron matter density in the sun has the values from lO^^^eV to 10~^^eV. Our 
result A ~ 10~^°eV, at which the four- neutrino resonance occurs, is very large. Therefore, 
one may not observe the sterile neutrino resonance realistically, even if it exists. To detect 
the fourth neutrino resonance, one needs matter of which density is about 10^^°eV. 

The average of the electron matter density in the earth is about 4.92 x lO^^^eV (see 
Appendix A). This is the same order as the matter density of the sun. So it will be difficult 
to find the resonance of the sterile neutrino near the earth. 

Figure 6 shows the transition probability Pss, which is the sterile neutrino surviving 
probability, as a function of the parameter tj for the (3-|-l)-scheme and the (2H-2)-scheme. 
From Fig. 6, the sterile neutrino transition occurs for 77 > 10~^. 



15 



The transition from an active neutrino to a sterile one may be observed if neutrino 
passes through the matter which has very high density. In the same condition, the reverse 
transition may occur. The above argument about the transition probability for four-neutrino 
oscillations tells that we may not realistically see the sterile neutrino resonance, even if the 
sterile neutrino exists. 
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APPENDIX A: ELECTRON NUMBER DENSITY IN THE NATURAL SYSTEM 

OF UNITS 

In this paper, we refer to the Ref. [16] for the physical constants. The electron number 
density N^, is connected with the matter density A by (22). Using the Fermi weak coupling 
constant Gp/ihcf = 1.17 x IQ-^fGeV-^], 

A[eY] = v^GpiVe = 1.27 x 10-^^[eV • m^] • Ne[l/m% (Al) 

where he = 197[MeV fm]. 

We discuss the matter density in the sun, using a solar model [17]. From the model, the 
electron number density depends on the distance R from the center of the sun : 

N^{x) = 98.19 Ar^e-^°-^5^ [cm-^], x = (A2) 

where Na = 6.02 x 10^^ and Rsun are the Avogadro constant and the radius of the sun, 
respectively. 

In Table I, the electron matter densities of the sun and the earth are listed. And we 
assume that the ratio of the electron mass to the nucleon mass is 1 : 2000, where the mass 
of electron is 9.11 x 10~^^[kg]. From Table I, the electron number density in the sun has the 
values from 10-^^eV to lO-^^eV. 
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TABLES 





A[eV] 


Ne [l/cm^' 


p [g/cm^] 


R/ Rsun — 

i?/i?su„ = 0.41 

R/ Rsun = 1 


1.06 xlO-^^ 

1.40 xlO"^''^ 
2.78 X 10-16 


5.91 xlO^^ 

7.82 xl023 
1.55 xlO^i 


108 

1.42 
2.82 xlO-2 


the earth 


5.44 X 10-13 


3.03 xl024 


5.52 



TABLE I. The electron matter densities of the sun and the earth. The value of A at 
X = R/Rsun = 0.41 corresponds to the average density in the sun. The mass density p is shown as 
p = 2000meAe, where rrie is the electron mass, and we put the ratio of the electron mass to the 
nucleon mass is 1 : 2000. We assume that the density in the earth is 5.52g/cm3. 
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FIGURES 
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FIG. 1. Two of mass patterns of four neutrino schemes. (a):(3+l)-scheme, (b):(2+2)-scheme. 




FIG. 2. The difference I Aa — Af, I (a, 6 = 1, 2, 3, 4, a / 6) as a function of the matter density A for 
a (3+l)-scheme, where 612 = f , 6*23 = f , 6*13 = 0, On = e, 624 = S = 0, 634 = 0, Am^i ~ IQ-^eV^, 
Am§2 =^ 10-3eV2, Am^i ~ leV^, S = 10 GeV and e = 0.1. 
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FIG. 3. The transition probabiHties Pee, Pes and Pss as a function of the matter density A for 
a (3+l)-scheme, where 612 = j, 623 = j, 6*13 = 0, 6*14 = e, 024 = S = 0, 6*34 = 0, Am|i ~ IQ-^eV^, 
Aml^ ~ IQ-^eV^ Amj^ ~ leV^, E = 10 GeV and L/E = 6.46 x lO^eV"^. The sohd and broken 
lines show the transition probabilities for e = 0.1,0, respectively. 
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FIG. 4. The differences |Aa — A^l (a, 6 = 1,2,3,4) as a function of the matter density A for 
a (2+2)-scheme, where = f , ^23 = f , ^12 = e, 6*13 = 0, ^24 = 0, ^34 = 0, Amji ~ leV^, 
Am|2 ^ 10-3eV2, Am|i ~ IQ-^eV^, E = 10 GeV and e = 0.1. 
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FIG. 5. The transition probabihties Pee, Pen and Pgs as a function of the matter density A 
for a (2+2)-schcmc, where 6^14 = f, 6*23 = f, 6*12 = e, 6*13 = 0, 6*24 = 0, 6*34 = 0, Amj^ ~ IcV^, 
Am|2 ^ lO-^eV^ Aml^ ~ IQ-^eV^, £; = 10 GeV and L/E = 6.46 x lO^eV-^. The solid and 
broken Unes show the transition probabihties for e = 0.1, 0, respectively. 
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P (3+1 scheme) P (2+2 scheme) 




FIG. 6. The surviving probabilities Pgg of the sterile neutrino transition as a function of the rate 

rj for a (3+l)-scheme and a (2+2)-scheme, where L/E = rj x 6.46 x lO^eV"^. We compare Pgs at 
two kinds of the matter densities for A = 1.40 x 10~^^[eV](the broken line) and A = 10~^''[eV](the 
solid line). 
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